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Time reversal mirrors work perfectly only for lossless wave propagation. Here, the performance 
of time-reversal mirrors is quantitatively defined, and the adverse effect of dissipation on their 
performance is investigated. The technique of exponential amplification is proposed to overcome 
the effect of dissipation in the case of uniform loss distributions, and, to some extent, in the case 
of non-uniform loss distributions. A numerical model of a star graph was employed to test the 
applicability of this technique on realizations with various random spatial distributions of loss. The 
numerical results are also verified by an experimental result from an electromagnetic time-reversal 
mirror. 



I. INTRODUCTION 

Time reversal (TR) mirrors can, under ideal circum- 
stances, precisely reconstruct a wave disturbance which 
happened at an earlier time, at any given later time. 
They have found numerous applications since the ear- 
liest experimental demonstrations of TR using acoustic 
waves pp. More recently, TR mirrors have been real- 
ized using electromagnetic waves [2j |3] expanding their 
range of applicability. They have applications in com- 
munications [4, 5 , imaging [6j[7], source localization [8], 
non-destructive evaluation [9j [TO] , and sensing [3j [TTJ [12] . 

TR mirrors are used to focus waves both in space and 
time. An ideal TR mirror consists of a wave source lo- 
cated inside a lossless medium that is completely enclosed 
by a surface of transceivers. A TR mirror operates in 
two steps. In step one, the transceivers record and ab- 
sorb the signal broadcast by the source. In step two, 
the transceivers rebroadcast time reversed versions of the 
recorded waves. The waves broadcast in step two even- 
tually focus on the location of the source and reconstruct 
a time reversed version of the original signal, which was 
broadcast in step one. This is possible because of the TR 
invariance property of the lossless wave equation. This 
ideal TR mirror uses the so called closed TR cavity [13] . 
However, it is not generally practical, for example, to 
build a closed TR cavity, whose interior is completely 
covered with transceivers. Practical TR mirrors have 
several limitations, which result in an imperfect recon- 
struction of the original signal. These limitations include 
i) limited spatial coverage by the transceivers, and ii) dis- 
sipation during the wave propagation (which breaks TR 
invariance of the wave equation) [TTJ [14] . 

The first limitation of TR mirrors can be overcome 
by the use of a reflecting wave chaotic cavity with par- 
tial spatial coverage and transceivers that have a long 
recording time [TTJ [14]. However, the limitation due to 
dissipation persists, and leads to increasing loss of infor- 
mation as the recording time increases. 

The technique of iterative TR was proposed to improve 
the performance of TR mirrors which suffer from the lim- 
itations of dissipation and incomplete spatial coverage of 



transceivers [T5] , However, the iterative technique may 
not converge fast enough, despite a recent result showing 
the tunability of its convergence [16] . A faster alternative 
to the iterative technique is the inverse filter technique 
[l5j[T7]. But, the inverse filter technique is computation- 
ally costly and potentially unstable [T5] . 

On the other hand, the technique of exponential am- 
plification was used to overcome the adverse effect of dis- 
sipation on TR based sensing techniques [TTJ [12] . The 
spatial range of these techniques is limited due to dis- 
sipation. The exponential amplification technique was 
used to extend the spatial range of the TR based sensing 
techniques. 

In this paper, another application of the exponential 
amplification technique is proposed. Unlike the work 
in Ref.[12], the application involves improving the per- 
formance of TR mirrors which suffer from dissipation. 
The performance of TR mirrors is defined as the nor- 
malized correlation between the original pulse and the 
time reversed version of the reconstructed pulse (after 
some post-processing). This paper quantitatively demon- 
strates that the exponential amplification technique mit- 
igates the adverse effects of dissipation on the perfor- 
mance of TR mirrors, for different spatial loss distribu- 
tions. The technique of exponential amplification can 
compete with the iterative TR technique in improving 
this performance of TR mirrors, to some extent. The ad- 
vantage of the exponential amplification technique is that 
it is faster to implement, as it does not rely on iterative 
steps. It is also cheaper computationally. 

In Sec. [H] the theory supporting the exponential am- 
plification technique, first discussed in Ref.[12 , is sum- 
marized. The performance of TR mirrors, which was 
not discussed in Ref.[12 , is also defined in Sec. [Hj In 
Sec. Ill the exponential amplification technique is used 
to improve the newly defined performance of an electro- 
magnetic TR mirror. In particular, it is experimentally 
shown that the technique can improve the performance of 
a TR mirror for the case of spatially uniform distributed 
loss. Sec. [IV] investigates the limits of applicability of 



the exponential amplification in systems with inhomo- 
geneous loss distributions. Sections [V] and VI present 
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further discussion and conclusions. 



II. THEORY 

The theory supporting the exponential amplification 
technique was originally presented in Ref.[12 j. In this 
section, the theory is summarized and the performance 
of TR mirrors is defined. To summarize the theory, a 
cavity with transfer function s(u) is considered. When 
a brief pulse, a(t), is broadcast into the cavity, the re- 
sponse is b(t). At a later time, a time reversed version of 
b(t) denoted b(—t) is broadcast back into the cavity pro- 
ducing reconstructed pulse denoted c{— t) at the original 
broadcasting antenna or transducer. 

Ref.[l2] provides a derivation showing that, if the cav- 
ity is lossless, a(t) is the same as a time reversed version 
of c(— t) denoted here c(t). Thus, the TR mirror works 
perfectly in a lossless cavity. 

Even if the cavity is lossy, it is shown that the TR 
mirror can still recover the original pulse in the special 
case of uniform loss. Suppose that the uniform loss is 
represented by a 1/e amplitude decay time r. In this 
case, the signal b(t) is exponentially amplified by e 2 ^ r , 
time reversed, and broadcast back into the cavity. b(—t) 
denotes the signal obtained after exponentially amplify- 
ing b(i), and time reversing it. Similarly, c(—t) denotes 
the reconstructed pulse obtained when b(—t) is broad- 
cast back into the cavity. Ref.[12] analytically predicts 
that c(i)e~ 2t l T is equal to a(t) for a uniform loss that is 
characterized by r. 

Here, we define a new signal c(t) which is equal to 
c(t)e~ 2t / T ". c(t) is called the corrected reconstructed 
pulse. In this paper, the performance (77) of the TR 
mirror is defined based on the normalized correlation be- 
tween the corrected reconstructed pulse (c(t)), and the 
original pulse (a(t)). 

,- , Eg "WW (1 ) 

Each of the signals a(t) and c(t) have a duration of T. In 
principle, 77 can range from —1 to 1. However, since the 
two signals are aligned to maximize 77, 77 ranges from 
(poor performance) to 1 (perfect reconstruction). Based 
on the result from Ref . [T2] . a(t) = c(t) if the system is 
uniformly lossy with a 1/e voltage decay time, r; there- 
fore, the performance of the TR mirror is expected to 
be perfect (i.e. 77 = 1). However, the analysis does not 
model the effect of noise and limited dynamic range on 
the exponential amplification technique. 

If the exponential amplification is not applied to &(£), 
the performance, 77, can still be calculated using the sig- 
nals a(t) and c(t), instead of using a(t) and c(t). But, if 
the system is lossy, and exponential amplification is not 
applied, then it is not expected that 77 = 1. 

In Sec. |III[ this theory, which assumes uniform loss dis- 
tribution, is successfully tested in an experimental system 



that approximates the case of uniform loss distribution. 
However, it is also expected that the technique might 
be useful in the case of a moderately non-uniform loss 
distribution. The applicability of the exponential ampli- 
fication technique to mitigate the effect of non-uniform 
loss is studied in Sec. [iVl 



III. OVERCOMING THE EFFECT OF 
SPATIALLY UNIFORM LOSS: EXPERIMENTAL 
TEST 

The electromagnetic TR mirror operates as shown in 
Fig. [l] The TR mirror involves an empty one cubic me- 
ter aluminum box that has scatterers inside it, and two 
ports in its walls. The aluminum walls of the cavity are 
constructed of a single material, and can be assumed to 
have a uniform spatial distribution of loss. The original 
pulse was broadcast into the cavity using port 1, and the 
response signal, which we call "sona", was collected us- 
ing port 2. The original pulse had a carrier frequency 
of 1GHz and a Gaussian envelope with a standard de- 
viation of Ins. The sona was recorded for about 6.5/is 
with a signal to noise ratio (SNR) greater than 1. Dur- 
ing the second step of the TR mirror operation, the sona 
signal was time reversed and broadcast back into the cav- 
ity using port 1; this made use of the spatial reciprocity 
property of the wave equation. Finally, a reconstructed 
pulse, which approximates a time reversed version of the 
original pulse was collected at port 2. The reconstructed 
pulse (and its temporal sidelobes) was collected over a 
lOfis duration. The first 6.5/is included a direct recording 
of the time reversed sona being injected into the cavity 
from port 1. The last 3.5 jis was a recording of the post- 
reconstructed pulse emerging from port 2, reverberating 
throughout the cavity. 

Based on Eq. [TJ 77 can be calculated as the normalized 
correlation between the following two aligned signals: i) 
the time reversed version of the reconstructed pulse with 
its temporal sidelobes (which was a 10/as long recording), 
and ii) the original pulse broadcast (by considering the 
Ins long Gaussian envelope to be located at t = 3.5us of 
a 10/is long signal of mostly zero voltage values). Note 
that these two signals are equivalent to the aligned a(t) 
and c(t) signals introduced in Sec. [Hj If the TR mirror 
were perfect, 77 would be 1. However, as can be seen 
in part 4 of Fig. [I] the reconstructed pulse has temporal 
sidelobes which result in 77 < 1. The sidelobes arise due to 
the limitations of practical TR mirrors discussed in Sec.[T| 
which include limited spatial coverage of transceivers and 
dissipation. 

Generally, the sona signal shows an exponential decay 
which is caused by dissipation and coupling losses. The 
experimental system discussed here is strongly under- 
coupled [18], which means that the effect of dissipation 
dominates the effect of coupling. Therefore, the 1/e volt- 
age decay constant (r) which is related to the unloaded 
quality factor can be readily determined by a line fit to 
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FIG. 1: Schematic of the electromagnetic time reversal (TR) 
mirror experiment, without exponential amplification. Dur- 
ing the first step of the TR mirror, the original pulse is broad- 
cast through antenna 1 (as shown in part 1), and the resulting 
sona is collected at antenna 2 (as shown in part 2). During 
the next step of the TR mirror, the time reversed sona is in- 
jected into the system at antenna 1 (as shown in part 3) to 
retrieve the reconstructed time reversed pulse at antenna 2 
using spatial reciprocity (as shown in part 4). Experimental 
data are shown for each step. 



FIG. 2: The reconstructed TR pulse obtained from a sona 
that was exponentially amplified with F — 2 (blue) had signif- 
icant sidelobes. The corrected reconstructed pulse (red) was 
obtained by multiplying the time reversed version of the re- 
constructed TR pulse by A(t — 3.5/is, F = —2). The corrected 
reconstructed pulse is displayed here after time reversing it. 
A(t,F) = 1 if t < 0. The inset shows a close up of the re- 
constructed pulses, which are the same before and after the 
correction. 



the log of the sona energy as a function of time. The ex- 
ponential amplification can then be applied to the sona 
signal before it is time reversed and broadcast into the 
cavity. The exponential amplification is carried out by 
multiplying the sona signal by the time dependent am- 
plifying function 

Ft 

A(t,F) = exp(—) (2) 
r 

where t is time in seconds, and F is an adjustable fac- 
tor. Assuming a precise determination of r, using F = 1 
compensates only for the effect of dissipation on the sona 
during step 1 of the TR mirror operation; whereas, using 
F = 2 compensates for the effect of dissipation on the 
sona during both step 1 and step 2 of the TR mirror op- 
eration [12]. The theoretical discussion in Sec. |TT| predicts 
that F = 2 to maximize rj. The exponential amplification 
is applied to the part of the sona whose signal to noise 
ratio is, at the very least, greater than 1 (i.e. the 6.5/is 
long sona in this experiment). In addition, the function 
A(t,F) is typically terminated by a smooth ramp-down 
function whose time span is at least as wide as the time 
duration of the original pulse, which is Ins; this prevents 
additional frequency components from entering into the 
sona. 

Suppose that the sona is exponentially amplified us- 
ing A(t, F), before it is time reversed and broadcast into 



the electromagnetic cavity. The reconstructed pulse ob- 
tained with F > will have more sidelobes than the 
reconstructed pulse obtained with F = 0. However, 
these sidelobes can be numerically suppressed by manip- 
ulating the reconstructed waveform as follows. Fig. [2] 
shows the reconstructed waveform obtained using sona 
that is amplified with A(t,F = 2) (in blue); the re- 
constructed waveform has unwanted temporal sidelobes 
which are a result of the exponential amplification ap- 
plied to the sona. These sidelobes were corrected as fol- 
lows. The 10/us long reconstructed waveform was time 
reversed so that the Ins long Gaussian pulse is located 
at about t = 3.5/is; this waveform was then multiplied by 
A(t — 3.5/is, F = —2) to get the corrected reconstructed 
waveform (i.e. c(t)) introduced in Sec. [ilj A(t,F) = 1 
for t < by convention. In Fig. [2j the corrected recon- 
structed waveform is time reversed and plotted (in red) 
to highlight its relationship with the reconstructed wave- 
form. The corrected reconstructed waveform is expected 
to have a better normalized correlation with the origi- 
nal pulse because of the exponential amplification with 
F = 2, as predicted in Sec. [TTJ 

Several values of F were used to carry out the op- 
eration of the TR mirror assisted by exponential am- 
plification. For each F value, the performance of the 
electromagnetic TR mirror, 77, is shown in Fig. [3j The 
performance of the TR mirror is enhanced the most by 
using exponential amplification with F « 2. In partic- 
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FIG. 3: The performance (77) of the experimental electro- 
magnetic TR mirror as a function of F parameter used to 
exponentially amplify the sona signal. The optimum value of 
the parameter is F — 2. 




FIG. 4: Schematic of the star graph model shows N transmis- 
sion lines and a driving transmission line, which are all con- 
nected in parallel. Each of the N lines (labeled by i) can have 
unique length (Li), characteristic admittance (Yd), frequency 
dependent propagation constant (7^(0;)), and termination re- 
flection coefficient (Ti). The driving line has a characteristic 
admittance of Y c d and zero length. 



ular, r] increases from « 63% at F = to « 73% at 
F = 2. This result agrees with the theoretical prediction 
in Sec. [TTJ The fact that the rj < 1 even for F « 2 can be 
explained by the fact that the TR mirror has only one 
recording channel, whereas the theory assumes that all 
the waves are captured. Thus, the exponential amplifica- 
tion technique only overcomes the limitation of the TR 
mirror associated with dissipation, and not the problem 
of insufficient spatial coverage of recorders. Noise, which 
is not considered in the theory, also affects the perfor- 
mance of TR mirrors as seen in Ref . [16] . 

As shown in Fig. [3j 77 « 63% if the exponential ampli- 
fication technique is not applied. The exponential am- 
plification with the optimum F = 2 parameter achieves 
r] ~ 73% for the experimental set up described. This 
compares with 77 « 78% which is achieved by the tun- 
able iterative technique on the same experimental set up 
[16]. Therefore, the iterative technique performs bet- 
ter than the exponential amplification technique when 
both of them use their respective optimum parameters. 
Nonetheless, the exponential amplification has an advan- 
tage because of its speed, and computational simplicity. 

The result in this section proves that the exponential 
amplification improves the performance, 77, of TR mirrors 
with uniform spatial loss. 



IV. OVERCOMING THE EFFECT OF 
SPATIALLY NON-UNIFORM LOSS: 
NUMERICAL TEST 

The exponential amplification technique is derived as- 
suming uniform loss distribution. The case of non- 



uniform loss distribution is better handled with tech- 
niques such as the inverse filter, and the iterative TR 
[T5] ITT] , which do not assume uniformity of loss. How- 
ever, it is expected that the exponential amplification can 
be applied in cavities with moderately inhomogeneous 
loss distributions. The applicability of the exponential 
amplification technique was studied as the inhomogene- 
ity of the spatial distribution of loss was increased. The 
simulations in this section were designed to supplement 



the experimental results in Sec III 



A. Simulation Setup 

The numerical model is an implementation of a star 
graph p~6j[T9]. The star graph is a quasi- ID chaotic sys- 
tem. It is modeled using a driving transmission line that 
feeds a number of other transmission lines, which are all 
connected to each other in parallel (see Fig. [4J. Each 
transmission line is terminated with a reflection load. 
The graph is a one port system. Thus, the original pulse 
is injected through the driving transmission line, and the 
sona is also collected from the same line. The original 
pulse has the same characteristics as the one used in the 



experiment described in Sec. Ill 



There are 500 transmission lines in the star graph set 
up, in addition to the driving line. Their length is given 
by Li = m for i ranging from 1 to 500. The driving 
line has a length of m. The characteristic admittances 
of the 500 lines are Y ci = IS. The driving line has a 
characteristic admittance of Y c d = X^=i°° — 5005; 
this set up eliminates prompt reflection of signals injected 
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through the driving line. The terminal reflection coeffi- 
cients, of the 500 lines are all set to 1. The dissipation 
is introduced through the frequency dependent propaga- 
tion constant of the lines, ji(uj) = i^+oli, where c is the 
speed of light and is the loss constant of line i. 

Oii specifies that the voltage waves decay on line i as 
e~ aiZ , where z is distance measured along line i in meters. 
The spatial inhomogeneity of loss on the star graph is 
modeled as follows. The of each of the 500 lines is 
randomly chosen from the probability density function 
(pdf), 

The pdf has two parameters, n and A, which define the 
mean (fi a = A + An) and standard deviation (a a = 
\J A(A + An)) of a values. This particular pdf is cho- 
sen for the following two reasons. First, the coefficient of 
variation (R = a a / u a ) can be easily varied while keeping 
ji a constant. Keeping the average spatial loss constant 
simplifies the problem, and helps us focus on the effect of 
increasing spatial loss inhomogeneity. The main motiva- 
tion for using the exponential amplification in the case of 
non-uniform loss is this: if the 1/e decay constants (i.e. 
r) of the different modes are not extremely different, we 
can use a single average r value. Thus, it is not interest- 
ing to vary the average loss (effectively r) here. In this 
set up ja a is chosen as (cr) _1 , where r = 1.5/is. The sec- 
ond reason to use this particular pdf is that its support is 
the set of positive numbers, which should be the case as 
a should always be positive on the passive transmission 
lines. Fig. [5] shows plots of the pdf for different values of 
R, where fi a w 0.002. 

The numerical values of the parameters of the star 
graph were chosen to achieve the following objective. The 
experimental result in Sec. |III| was based on an under- 
coupled cavity. An under-coupled cavity is characterized 
by the domination of energy loss due to dissipation over 
energy loss due to coupling [15] . If the cavity is not 
under-coupled, the advantage of the exponential amplifi- 
cation technique cannot be easily seen because there are 
no strong internal dissipation effects to be compensated 
in an over-coupled cavity. To accomplish under-coupling, 
500 lines were used in the star graph to increase the back- 
reflection coefficient of the waves that are trying to leave 
the star graph. This forced the waves to reverberate in- 
side the star graph longer, which decreased the coupling 
loss by two orders of magnitude compared to the dissi- 
pation loss. 
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FIG. 5: The probability density function (pdf) in Eq. [3] is 
plotted for different R = cr a /fi a values where \i a ~ 0.002. R 
is a measure of loss inhomogeneity. The pdf is plotted for 
R = 0.2 (blue), R = 0.5 (black), and R = 1 (red). 
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FIG. 6: The performance (i.e. 77) of the TR mirror in a star 
graph model as a function of F parameter for various degrees 
of loss inhomogeneities characterized by R. Loss inhomogene- 
ity increases with R. The optimum F value is always around 
F — 2 for the r value that is well determined from the sona. 



B. Simulation Results 

A TR mirror was implemented on the star graph by 
broadcasting the original pulse, collecting a sona, expo- 
nentially amplifying the sona with A(t,F), time revers- 
ing the amplified sona, and broadcasting it back into the 
star graph. The reconstructed pulse was collected, and 



it was time reversed and multiplied by A(t — t p ,—F). 
Here, t p is the time when the Ins Gaussian pulse was 
expected within the time reversed version of the recon- 
structed waveform. The signals were collected over a 
span of lOus without noise. 

The performance (77) of the TR mirror is plotted versus 
F as shown in Fig. [6] This entire process was repeated 
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FIG. 7: The performance of the TR mirror in a star graph 
model is plotted as a function of the loss inhomogeneity, R, for 
the cases of: i) no exponential amplification (shown in black) , 
and ii) optimal exponential amplification, which is always F — 
2 (shown in red). It is clear that exponential amplification is 
advantageous for uniform and moderately non-uniform spatial 
loss distributions. 



for star graphs with different degree of loss inhomogene- 
ity characterized by R. R = represents uniform spatial 
loss distribution. As R increases, the loss inhomogene- 
ity increases. For each R value, 25 realizations of the 
star graph were generated; hence, there are error bars 
included on the 77 vs F plots for each R value, as shown 
in Fig. [6j 

For each R value shown in Fig. [6j the maximum 77 was 
taken over all F values attempted; it turns out that the 
maximum 77 is always achieved for F — 2. The maxi- 
mum 77 achieved is plotted as a function of the loss inho- 
mogeneity (R) in Fig. [7| (shown in red). The 77 achieved 
without using exponential amplification (i.e. F = 0) is 
also plotted as a function of R in Fig. [7] (shown in black). 
Fig. [7] demonstrates that the exponential amplification 
technique improves 77 significantly, for small values of R, 
which are close to uniform loss distributions. The ex- 
ponential amplification is still applicable in cavities with 
more inhomogeneous loss distributions. However, its ef- 
fectiveness declines with increasing loss inhomogeneity as 
expected. 

From Fig. [7| we can see that 77 is not affected signif- 
icantly by R if exponential amplification is not applied. 
Here, it is worth repeating the fact that the average spa- 
tial loss in the star graph is constant for all the R values 
considered. On the other hand, the fact that 77 is not 
restored to 1 even for R = and F — 2 can be explained 
by the fact that the star graph is designed to be under- 
coupled, to match the nature of the experimental cavity 
in Sec. |III| The waves in the star graph will continue to 



reverberate for longer than the 10/is sona recording time. 
Thus, the problem is not just dissipation which can be 
compensated by the exponential amplification technique, 
but the problem here is also the lack of sufficient tempo- 
ral coverage, which has previously been shown to affect 
the quality of TR focusing [20] . 



V. DISCUSSION 

The tunable iterative time-reversal technique over- 
comes the adverse effect of dissipation on the real-time 
spatiotemporal wave focusing of time reversal mirrors 
[T6] . The focusing was quantified using the ratio of pulse 
to temporal sidelobe energy. In this paper, we are less in- 
terested in achieving a real-time spatiotemporal wave fo- 
cusing using time reversal mirrors. Instead we are mainly 
interested in digitizing the reconstructed time reversed 
pulse, and numerically processing it to suppress the tem- 
poral sidelobes of the pulse. Suppressing the tempo- 
ral sidelobes improves the fidelity of the reconstructed 
pulse with the time reversed version of the original pulse 
(i.e. a(—t)). The temporal sidelobes of the reconstructed 
pulse can be significant, and they are created as a byprod- 
uct of the exponential amplification of the sona. Fig. [2] 
shows a digitized reconstructed pulse before and after its 
temporal sidelobes are numerically suppressed. The ca- 
pability to systematically suppress these temporal side- 
lobes after the digitization of the reconstructed pulse can 
be useful. For instance, if the shape of the reconstructed 
pulse is encoding information in a communication appli- 
cation, it is important to be able to improve the fidelity 
of the reconstructed pulse (which is measured by 77) by 
suppressing its temporal sidelobes. 



VI. CONCLUSION 

Exponential amplification improves the performance of 
time reversal mirrors best if the loss is uniformly dis- 
tributed in space. However, even under conditions in 
which the loss is not uniformly distributed, the expo- 
nential amplification mitigates the adverse effect of non- 
uniform loss on the performance of time reversal mirrors. 
The exponential amplification works best when the theo- 
retically predicted optimum parameter of F = 2 is used. 
This technique may be advantageous to use in time re- 
versal mirror based communication applications because 
even if the reconstructed pulse will have significant tem- 
poral sidelobes during the recording, the processed re- 
constructed pulse will closely replicate the shape of the 
original pulse. 

This work is supported by ONR MURI grant 
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